Analysis of polygenic selection in purebred and crossbred pig genomes using generation proxy selection mapping.

BACKGROUND: Artificial selection on quantitative traits using breeding values and selection indices in commercial livestock breeding populations causes changes in allele frequency over time at hundreds or thousands of causal loci and the surrounding genomic regions. In population genetics, this type of selection is called polygenic selection. Researchers and managers of pig breeding programs are motivated to understand the genetic basis of phenotypic diversity across genetic lines, breeds, and populations using selection mapping analyses. Here, we applied generation proxy selection mapping (GPSM), a genome-wide association analysis of single nucleotide polymorphism (SNP) genotypes (38,294-46,458 markers) of birth date, in four pig populations (15,457, 15,772, 16,595 and 8447 pigs per population) to identify loci responding to artificial selection over a period of five to ten years. Gene-drop simulation analyses were conducted to provide context for the GPSM results. Selected loci within and across each population of pigs were compared in the context of swine breeding objectives. RESULTS: The GPSM identified 49 to 854 loci as under selection (Q-values less than 0.10) across 15 subsets of pigs based on combinations of populations. The number of significant associations increased when data were pooled across populations. In addition, several significant associations were identified in more than one population. These results indicate concurrent selection objectives, similar genetic architectures, and shared causal variants responding to selection across these pig populations. Negligible error rates (less than or equal to 0.02%) of false-positive associations were found when testing GPSM on gene-drop simulated genotypes, suggesting that GPSM distinguishes selection from random genetic drift in actual pig populations. CONCLUSIONS: This work confirms the efficacy and the negligible error rates of the GPSM method in detecting selected loci in commercial pig populations. Our results suggest shared selection objectives and genetic architectures across swine populations. The identified polygenic selection highlights loci that are important to swine production.

Genotyping and phenotyping strategies for genetic improvement of meat quality and carcass composition in swine.

BACKGROUND: Meat quality and composition traits have become valuable in modern pork production; however, genetic improvement has been slow due to high phenotyping costs. Combining genomic information with multi-trait indirect selection based on cheaper indicator traits is an alternative for continued cost-effective genetic improvement. METHODS: Data from an ongoing breeding program were used in this study. Phenotypic and genomic information was collected on three-way crossbred and purebred Duroc animals belonging to 28 half-sib families. We applied different methods to assess the value of using purebred and crossbred information (both genomic and phenotypic) to predict expensive-to-record traits measured on crossbred individuals. Estimation of multi-trait variance components set the basis for comparing the different scenarios, together with a fourfold cross-validation approach to validate the phenotyping schemes under four genotyping strategies. RESULTS: The benefit of including genomic information for multi-trait prediction depended on the breeding goal trait, the indicator traits included, and the source of genomic information. While some traits benefitted significantly from genotyping crossbreds (e.g., loin intramuscular fat content, backfat depth, and belly weight), multi-trait prediction was advantageous for some traits even in the absence of genomic information (e.g., loin muscle weight, subjective color, and subjective firmness). CONCLUSIONS: Our results show the value of using different sources of phenotypic and genomic information. For most of the traits studied, including crossbred genomic information was more beneficial than performing multi-trait prediction. Thus, we recommend including crossbred individuals in the reference population when these are phenotyped for the breeding objective.

Imputation of non-genotyped F1 dams to improve genetic gain in swine crossbreeding programs.

This study investigated using imputed genotypes from non-genotyped animals which were not in the pedigree for the purpose of genetic selection and improving genetic gain for economically relevant traits. Simulations were used to mimic a 3-breed crossbreeding system that resembled a modern swine breeding scheme. The simulation consisted of three purebred (PB) breeds A, B, and C each with 25 and 425 mating males and females, respectively. Males from A and females from B were crossed to produce AB females (n = 1,000), which were crossed with males from C to produce crossbreds (CB; n = 10,000). The genome consisted of three chromosomes with 300 quantitative trait loci and ~9,000 markers. Lowly heritable reproductive traits were simulated for A, B, and AB (h2 = 0.2, 0.2, and 0.15, respectively), whereas a moderately heritable carcass trait was simulated for C (h2 = 0.4). Genetic correlations between reproductive traits in A, B, and AB were moderate (rg = 0.65). The goal trait of the breeding program was AB performance. Selection was practiced for four generations where AB and CB animals were first produced in generations 1 and 2, respectively. Non-genotyped AB dams were imputed using FImpute beginning in generation 2. Genotypes of PB and CB were used for imputation. Imputation strategies differed by three factors: 1) AB progeny genotyped per generation (2, 3, 4, or 6), 2) known or unknown mates of AB dams, and 3) genotyping rate of females from breeds A and B (0% or 100%). PB selection candidates from A and B were selected using estimated breeding values for AB performance, whereas candidates from C were selected by phenotype. Response to selection using imputed genotypes of non-genotyped animals was then compared to the scenarios where true AB genotypes (trueGeno) or no AB genotypes/phenotypes (noGeno) were used in genetic evaluations. The simulation was replicated 20 times. The average increase in genotype concordance between unknown and known sire imputation strategies was 0.22. Genotype concordance increased as the number of genotyped CB increased with little additional gain beyond 9 progeny. When mates of AB were known and more than 4 progeny were genotyped per generation, the phenotypic response in AB did not differ (P > 0.05) from trueGeno yet was greater (P < 0.05) than noGeno. Imputed genotypes of non-genotyped animals can be used to increase performance when 4 or more progeny are genotyped and sire pedigrees of CB animals are known.

In swine breeding, phenotypic information is often gathered from elite purebred (PB) breeding stock and occasionally terminal crossbred animals (CB). Using economically relevant traits expressed by dams of CB (F1) in genetic evaluations is not common due to the lack of pedigree and/or genomic relationships to relate phenotypes of F1 to PB selection candidates. Since swine often have large litters, this study aimed to develop strategies to incorporate phenotypes of F1 into genetic evaluations by imputing F1 genotypes. Using simulation, we investigated the impact of CB pedigree completeness, the number of CB genotyped progeny, the number of parities (and thus mates) a F1 had, and genomic diversity in PB breeds on imputation accuracy and the response to selection in F1 performance. When mates of F1 were in the pedigree and 4 or more CB progeny were genotyped per generation, imputation accuracy was high and the phenotypic response in F1 did not differ compared to when true F1 genotypes were used. Our results show that imputed genotypes can be used to increase performance in swine breeding programs, but the magnitude depends upon the number of CB progeny genotyped, the number of F1 mates, and the completeness of the pedigree.

Effects of boar sperm antioxidant supplementation on fertility.

Use of cooled stored boar semen for artificial insemination accounts for 99% of inseminations in pork production sectors worldwide. The aim of the present study was to examine use of GameteGuard®-CP, a natural plant derived blend of antioxidants, supplemented in one commercially available semen extender on sperm quality, conception rate, farrowing rate, and litter size. Ejaculates from 16 commercial Duroc boars were used in a split-ejaculate single-sire artificial insemination study. Sperm were evaluated for motility, viability, acrosome integrity, membrane stability, mitochondrial membrane potential, DNA intactness, and total antioxidant reactivity on both day 1 and 4 subsequent to ejaculate collection. Sows were inseminated after a treatment regimen was imposed to synchronize time of ovulation (n = 1476) among females using sperm collected 4 days subsequent to ejaculate collection extended in AndroStar Plus® or GameteGuard®-CP supplemented extender. At day 4 post-ejaculate collection, there was maintenance of sperm viability, membrane stability, acrosome intactness, mitochondrial membrane potential, and DNA intactness (P < 0.01) when there was use of GameteGuard®-CP treatments. When control samples were used, there were lesser (P < 0.01) values from day 1-4 for all variables evaluated. Conception rate and farrowing rates were 7.4% and 9.7% greater, respectively (P < 0.05), but there was no difference as a result of GameteGuard®-CP-treatment in mean litter outcomes. The results in the present study indicate antioxidant supplementation using GameteGuard®-CP maintains sperm quality during cool sperm storage resulting in improved conception and farrowing rates when using semen stored in the cooled state for as long as 4 days.

Microbial composition differs between production systems and is associated with growth performance and carcass quality in pigs.

BACKGROUND: The role of the microbiome in livestock production has been highlighted in recent research. Currently, little is known about the microbiome's impact across different systems of production in swine, particularly between selection nucleus and commercial populations. In this paper, we investigated fecal microbial composition in nucleus versus commercial systems at different time points. RESULTS: We identified microbial OTUs associated with growth and carcass composition in each of the two populations, as well as the subset common to both. The two systems were represented by individuals with sizeable microbial diversity at weaning. At later times microbial composition varied between commercial and nucleus, with species of the genus Lactobacillus more prominent in the nucleus population. In the commercial populations, OTUs of the genera Lactobacillus and Peptococcus were associated with an increase in both growth rate and fatness. In the nucleus population, members of the genus Succinivibrio were negatively correlated with all growth and carcass traits, while OTUs of the genus Roseburia had a positive association with growth parameters. Lactobacillus and Peptococcus OTUs showed consistent effects for fat deposition and daily gain in both nucleus and commercial populations. Similarly, OTUs of the Blautia genus were positively associated with daily gain and fat deposition. In contrast, an increase in the abundance of the Bacteroides genus was negatively associated with growth performance parameters. CONCLUSIONS: The current study provides a first characterization of microbial communities' value throughout the pork production systems. It also provides information for incorporating microbial composition into the selection process in the quest for affordable and sustainable protein production in swine.

Gut microbiome mediates host genomic effects on phenotypes: a case study with fat deposition in pigs.

A large number of studies have highlighted the importance of gut microbiome composition in shaping fat deposition in mammals. Several studies have also highlighted how host genome controls the abundance of certain species that make up the gut microbiota. We propose a systematic approach to infer how the host genome can control the gut microbiome, which in turn contributes to the host phenotype determination. We implemented a mediation test that can be applied to measured and latent dependent variables to describe fat deposition in swine (Sus scrofa). In this study, we identify several host genomic features having a microbiome-mediated effects on fat deposition. This demonstrates how the host genome can affect the phenotypic trait by inducing a change in gut microbiome composition that leads to a change in the phenotype. Host genomic variants identified through our analysis are different than the ones detected in a traditional genome-wide association study. In addition, the use of latent dependent variables allows for the discovery of additional host genomic features that do not show a significant effect on the measured variables. Microbiome-mediated host genomic effects can help understand the genetic determination of fat deposition. Since their contribution to the overall genetic variance is usually not included in association studies, they can contribute to filling the missing heritability gap and provide further insights into the host genome - gut microbiome interplay. Further studies should focus on the portability of these effects to other populations as well as their preservation when pro-/pre-/anti-biotics are used (i.e. remediation).

Microbiability of meat quality and carcass composition traits in swine.

The impact of gut microbiome composition was investigated at different stages of production (weaning, Mid-test and Off-test) on meat quality and carcass composition traits of 1,123 three-way crossbred pigs. Data were analysed using linear mixed models which included the fixed effects of dam line, contemporary group and gender as well as the random effects of pen, animal and microbiome information at different stages. The contribution of the microbiome to all traits was prominent although it varied over time, increasing from weaning to Off-test for most traits. Microbiability estimates of carcass composition traits were greater than that of meat quality traits. Among all of the traits analysed, belly weight (BEL) had a higher microbiability estimate (0.29 ± 0.04). Adding microbiome information did not affect the estimates of genomic heritability of meat quality traits but affected the estimates of carcass composition traits. Fat depth had a greater decrease (10%) in genomic heritability at Off-test. High microbial correlations were found among different traits, particularly with traits related to fat deposition with a decrease in the genomic correlation up to 20% for loin weight and BEL. This suggested that genomic correlation was partially contributed by genetic similarity of microbiome composition. The results indicated that better understanding of microbial composition could aid the improvement of complex traits, particularly the carcass composition traits in swine by inclusion of microbiome information in the genetic evaluation process.

Genomics of Heat Tolerance in Reproductive Performance Investigated in Four Independent Maternal Lines of Pigs.

Improving swine climatic resilience through genomic selection has the potential to minimize welfare issues and increase the industry profitability. The main objective of this study was to investigate the genetic and genomic determinism of tolerance to heat stress in four independent purebred populations of swine. Three female reproductive traits were investigated: total number of piglets born (TNB), number of piglets born alive (NBA) and average birth weight (ABW). More than 80,000 phenotypic and 12,000 genotyped individuals were included in this study. Genomic random-regression models were fitted regressing the phenotypes of interest on a set of 95 environmental covariates extracted from public weather station records. The models yielded estimates of (genomic) reactions norms for individual pigs, as indicator of heat tolerance. Heat tolerance is a heritable trait, although the heritabilities are larger under comfortable than heat-stress conditions (larger than 0.05 vs. 0.02 for TNB; 0.10 vs. 0.05 for NBA; larger than 0.20 vs. 0.10 for ABW). TNB showed the lowest genetic correlation (-38%) between divergent climatic conditions, being the trait with the strongest impact of genotype by environment interaction, while NBA and ABW showed values slightly negative or equal to zero reporting a milder impact of the genotype by environment interaction. After estimating genetic parameters, a genome-wide association study was performed based on the single-step GBLUP method. Heat tolerance was observed to be a highly polygenic trait. Multiple and non-overlapping genomic regions were identified for each trait based on the genomic breeding values for reproductive performance under comfortable or heat stress conditions. Relevant regions were found on chromosomes (SSC) 1, 3, 5, 6, 9, 11, and 12, although there were important regions across all autosomal chromosomes. The genomic region located on SSC9 appears to be of particular interest since it was identified for two traits (TNB and NBA) and in two independent populations. Heat tolerance based on reproductive performance indicators is a heritable trait and genetic progress for heat tolerance can be achieved through genetic or genomic selection. Various genomic regions and candidate genes with important biological functions were identified, which will be of great value for future functional genomic studies.

Modeling host-microbiome interactions for the prediction of meat quality and carcass composition traits in swine.

BACKGROUND: The objectives of this study were to evaluate genomic and microbial predictions of phenotypes for meat quality and carcass traits in swine, and to evaluate the contribution of host-microbiome interactions to the prediction. Data were collected from Duroc-sired three-way crossbred individuals (n = 1123) that were genotyped with a 60 k SNP chip. Phenotypic information and fecal 16S rRNA microbial sequences at three stages of growth (Wean, Mid-test, and Off-test) were available for all these individuals. We used fourfold cross-validation with animals grouped based on sire relatedness. Five models with three sets of predictors (full, informatively reduced, and randomly reduced) were evaluated. 'Full' included information from all genetic markers and all operational taxonomic units (OTU), while 'informatively reduced' and 'randomly reduced' represented a reduced number of markers and OTU based on significance preselection and random sampling, respectively. The baseline model included the fixed effects of dam line, sex and contemporary group and the random effect of pen. The other four models were constructed by including only genomic information, only microbiome information, both genomic and microbiome information, and microbiome and genomic information and their interaction. RESULTS: Inclusion of microbiome information increased predictive ability of phenotype for most traits, in particular when microbiome information collected at a later growth stage was used. Inclusion of microbiome information resulted in higher accuracies and lower mean squared errors for fat-related traits (fat depth, belly weight, intramuscular fat and subjective marbling), objective color measures (Minolta a*, Minolta b* and Minolta L*) and carcass daily gain. Informative selection of markers increased predictive ability but decreasing the number of informatively reduced OTU did not improve model performance. The proportion of variation explained by the host-genome-by-microbiome interaction was highest for fat depth (~ 20% at Mid-test and Off-test) and shearing force (~ 20% consistently at Wean, Mid-test and Off-test), although the inclusion of the interaction term did not increase the accuracy of predictions significantly. CONCLUSIONS: This study provides novel insight on the use of microbiome information for the phenotypic prediction of meat quality and carcass traits in swine. Inclusion of microbiome information in the model improved predictive ability of phenotypes for fat deposition and color traits whereas including a genome-by-microbiome term did not improve prediction accuracy significantly.

Heritability and genome-wide association of swine gut microbiome features with growth and fatness parameters.

Despite recent efforts to characterize longitudinal variation in the swine gut microbiome, the extent to which a host's genome impacts the composition of its gut microbiome is not yet well understood in pigs. The objectives of this study were: i) to identify pig gut microbiome features associated with growth and fatness, ii) to estimate the heritability of those features, and, iii) to conduct a genome-wide association study exploring the relationship between those features and single nucleotide polymorphisms (SNP) in the pig genome. A total of 1,028 pigs were characterized. Animals were genotyped with the Illumina PorcineSNP60 Beadchip. Microbiome samples from fecal swabs were obtained at weaning (Wean), at mid-test during the growth trial (MidTest), and at the end of the growth trial (OffTest). Average daily gain was calculated from birth to week 14 of the growth trial, from weaning to week 14, from week 14 to week 22, and from week 14 to harvest. Backfat and loin depth were also measured at weeks 14 and 22. Heritability estimates (±SE) of Operational Taxonomic Units ranged from 0.025 (±0.0002) to 0.139 (±0.003), from 0.029 (±0.003) to 0.289 (±0.004), and from 0.025 (±0.003) to 0.545 (±0.034) at Wean, MidTest, and OffTest, respectively. Several SNP were significantly associated with taxa at the three time points. These SNP were located in genomic regions containing a total of 68 genes. This study provides new evidence linking gut microbiome composition with growth and carcass traits in swine, while also identifying putative host genetic markers associated with significant differences in the abundance of several prevalent microbiome features.

Genetic Parameters for Tolerance to Heat Stress in Crossbred Swine Carcass Traits.

Data for loin and backfat depth, as well as carcass growth of 126,051 three-way crossbred pigs raised between 2015 and 2019, were combined with climate records of air temperature, relative humidity, and temperature-humidity index. Environmental covariates with the largest impact on the studied traits were incorporated in a random regression model that also included genomic information. Genetic control of tolerance to heat stress and the presence of genotype by environment interaction were detected. Its magnitude was more substantial for loin depth and carcass growth, but all the traits studied showed a different impact of heat stress and different magnitude of genotype by environment interaction. For backfat depth, heritability was larger under comfortable conditions (no heat stress), as compared to heat stress conditions. Genetic correlations between extreme values of environmental conditions were lower (∼0.5 to negative) for growth and loin depth. Based on the solutions obtained from the model, sires were ranked on their breeding value for general performance and tolerance to heat stress. Antagonism between overall performance and tolerance to heat stress was moderate. Still, the models tested can provide valuable information to identify genetic material that is resilient and can perform equally when environmental conditions change. Overall, the results obtained from this study suggest the existence of genotype by environment interaction for carcass traits, as a possible genetic contributor to heat tolerance in swine.

Genome-wide association study for carcass quality traits and growth in purebred and crossbred pigs1.

Carcass quality traits such as back fat (BF), loin depth (LD), and ADG are of extreme economic importance for the swine industry. This study aimed to (i) estimate the genetic parameters for such traits and (ii) conduct a single-step genome-wide association study (ssGWAS) to identify genomic regions that affect carcass quality and growth traits in purebred (PB) and three-way crossbred (CB) pigs. A total of 28,497 PBs and 135,768 CBs pigs were phenotyped for BF, LD, and ADG. Of these, 4,857 and 3,532 were genotyped using the Illumina PorcineSNP60K Beadchip. After quality control, 36,328 SNPs were available and were used to perform an ssGWAS. A bootstrap analysis (n = 1,000) and a signal enrichment analysis were performed to declare SNP significance. Genome regions were based on the variance explained by significant 10-SNP sliding windows. Estimates of PB heritability (SE) were 0.42 (0.019) for BF, 0.39 (0.020) for LD, and 0.35 (0.021) for ADG. Estimates of CB heritability were 0.49 (0.042) for BF, 0.27 (0.029) for LD, and 0.12 (0.021) for ADG. Genetic correlations (SE) across the two populations were 0.81 (0.02), 0.79 (0.04), and 0.56 (0.05), for BF, LD, and ADG, respectively. The variance explained by significant regions for each trait in PBs ranged from 1.51% to 1.35% for BF, from 4.02% to 3.18% for LD, and from 2.26% to 1.45% for ADG. In CBs, the variance explained by significant regions ranged from 1.88% to 1.37% for BF, from 1.29% to 1.23% for LD, and from 1.54% to 1.32% for ADG. In this study, we have described regions of the genome that determine carcass quality and growth traits of PB and CB pigs. These results provide evidence that there are overlapping and nonoverlapping regions in the genome influencing carcass quality and growth traits in PBs and three-way CB pigs.

Genetic parameters of meat quality, carcass composition, and growth traits in commercial swine.

Swine industry breeding goals are mostly directed towards meat quality and carcass traits due to their high economic value. Yet, studies on meat quality and carcass traits including both phenotypic and genotypic information remain limited, particularly in commercial crossbred swine. The objectives of this study were to estimate the heritabilities for different carcass composition traits and meat quality traits and to estimate the genetic and phenotypic correlations between meat quality, carcass composition, and growth traits in 2 large commercial swine populations: The Maschhoffs LLC (TML) and Smithfield Premium Genetics (SPG), using genotypes and phenotypes data. The TML data set consists of 1,254 crossbred pigs genotyped with 60K SNP chip and phenotyped for meat quality, carcass composition, and growth traits. The SPG population included over 35,000 crossbred pigs phenotyped for meat quality, carcass composition, and growth traits. For TML data sets, the model included fixed effects of dam line, contemporary group (CG), gender, as well as random additive genetic effect and pen nested within CG. For the SPG data set, fixed effects included parity, gender, and CG, as well as random additive genetic effect and harvest group. Analyses were conducted using BLUPF90 suite of programs. Univariate and bivariate analyses were implemented to estimate heritabilities and correlations among traits. Primal yield traits were uniquely created in this study. Heritabilities [high posterior density interval] of meat quality traits ranged from 0.08 [0.03, 0.16] for pH and 0.08 [0.03, 0.1] for Minolta b* to 0.27 [0.22, 0.32] for marbling score, except intramuscular fat with the highest estimate of 0.52 [0.40, 0.62]. Heritabilities of primal yield traits were higher than that of primal weight traits and ranged from 0.17 [0.13, 0.25] for butt yield to 0.45 [0.36, 0.55] for ham yield. The genetic correlations of meat quality and carcass composition traits with growth traits ranged from moderate to high in both directions. High genetic correlations were observed for male and female for all traits except pH. The genetic parameter estimates of this study indicate that a multitrait approach should be considered for selection programs aimed at meat quality and carcass composition in commercial swine populations.

Predicting Growth and Carcass Traits in Swine Using Microbiome Data and Machine Learning Algorithms.

In this paper, we evaluated the power of microbiome measures taken at three time points over the growth test period (weaning, 15 and 22 weeks) to foretell growth and carcass traits in 1039 individuals of a line of crossbred pigs. We measured prediction accuracy as the correlation between actual and predicted phenotypes in a five-fold cross-validation setting. Phenotypic traits measured included live weight measures and carcass composition obtained during the trial as well as at slaughter. We employed a null model excluding microbiome information as a baseline to assess the increase in prediction accuracy stemming from the inclusion of operational taxonomic units (OTU) as predictors. We further contrasted performance of models from the Bayesian alphabet (Bayesian Lasso) as well machine learning approaches (Random Forest and Gradient Boosting) and semi-parametric kernel models (Reproducing Kernel Hilbert space). In most cases, prediction accuracy increased significantly with the inclusion of microbiome data. Accuracy was more substantial with the inclusion of microbiome information taken at weeks 15 and 22, with values ranging from approximately 0.30 for loin traits to more than 0.50 for back fat. Conversely, microbiome composition at weaning resulted in most cases in marginal gains of prediction accuracy, suggesting that later measures might be more useful to include in predictive models. Model choice affected predictions marginally with no clear winner for any model/trait/time point. We, therefore, suggest average prediction across models as a robust strategy in fitting microbiome information. In conclusion, microbiome composition can effectively be used as a predictor of growth and composition traits, particularly for fatness traits. The inclusion of OTU predictors could potentially be used to promote fast growth of individuals while limiting fat accumulation. Early microbiome measures might not be good predictors of growth and OTU information might be best collected at later life stages. Future research should focus on the inclusion of both microbiome as well as host genome information in predictions, as well as the interaction between the two. Furthermore, the influence of the microbiome on feed efficiency as well as carcass and meat quality should be investigated.

The effect of a porcine reproductive and respiratory syndrome outbreak on genetic parameters and reaction norms for reproductive performance in pigs1.

The objective of this study was to estimate genetic parameters of antibody response and reproductive traits after exposure to porcine reproductive and respiratory syndrome virus. Blood samples were taken approximately 60 d after the outbreak. Antibody levels were quantified as the sample-to-positive ratio (S/P ratio) using a fluorescent microsphere assay. Reproductive traits included total number born (TNB), number born alive (NBA), number stillborn (NSB), number mummified (NBM), and number born dead (NBD). Mortality traits were log transformed for genetic analyses. Data were split into prior, during, and after the disease outbreak phases using visual appraisal of the estimates of farm-year-week effects for each reproductive trait. For NBA, data from all phases were combined into a reaction norm analysis with regression on estimates of farm-year-week effects for NBA. Heritability for S/P ratio was estimated at 0.17 ± 0.05. Heritability estimates for reproduction traits were all low and were lower during the outbreak for NBA but greater for mortality traits. TNB was not greatly affected during the outbreak, as many sows that farrowed during the outbreak were mated prior to the outbreak. Heritability for TNB decreased from 0.13 (prior) to 0.08 (after). Genetic correlation estimates between prior to and during the outbreak were high for TNB (0.86 ± 0.23) and NBA (0.98 ± 0.38) but lower for mortality traits: 0.65 ± 0.43, -0.42 ± 0.55, and 0.29 ± 1.39 for LNSB, LNBM, and LNBD, respectively. TNB prior to and after the outbreak had a lower genetic correlation (0.32 ± 0.33). In general, genetic correlation estimates of S/P ratio with reproductive performance during the outbreak were below 0.20 in absolute value, except for LNSB (-0.73 ± 0.29). Based on the reaction norm model, estimates of genetic correlations between the intercept and slope terms ranged from 0.24 ± 0.50 to 0.54 ± 0.35 depending on the parameterization used, indicating that selection for the intercept may result in indirect selection for steeper slopes, and thus, less resilient animals. In general, estimates of genetic correlations between farm-year-week effect classes based on the reaction norm model resembled estimates of genetic correlations from the multivariate analysis. Overall, compared to previous studies, antibody S/P ratios showed a lower heritability (0.17 ± 0.05) and low genetic correlations with reproductive performance during a porcine reproductive and respiratory syndrome outbreak, except for the LNSB.

Host contributes to longitudinal diversity of fecal microbiota in swine selected for lean growth.

BACKGROUND: In pigs, gut bacteria have been shown to play important roles in nutritional, physiological, and immunological processes in the host. However, the contribution of their metagenomes or part of them, which are normally reflected by fragments of 16S rRNA-encoding genes, has yet to be fully investigated. RESULTS: Fecal samples, collected from a population of crossbred pigs at three time points, including weaning, week 15 post weaning (hereafter "week 15"), and end-of-feeding test (hereafter "off-test"), were used to evaluate changes in the composition of the fecal microbiome of each animal over time. This study used 1205, 1295, and 1283 samples collected at weaning, week 15, and off-test, respectively. There were 1039 animals that had samples collected at all three time points and also had phenotypic records on back fat thickness (BF) and average daily body weight gain (ADG). Firmicutes and Bacteroidetes were the most abundant phyla at all three time points. The most abundant genera at all three time points included Clostridium, Escherichia, Bacteroides, Prevotella, Ruminococcus, Fusobacterium, Campylobacter, Eubacterium, and Lactobacillus. Two enterotypes were identified at each time point. However, only enterotypes at week 15 and off-test were significantly associated with BF. We report herein two novel findings: (i) alpha diversity and operational taxonomic unit (OTU) richness were moderately heritable at week 15, h2 of 0.15 ± 0.06 to 0.16 ± 0.07 and 0.23 ± 0.09 to 0.26 ± 0.08, respectively, as well as at off-test, h2 of 0.20 ± 0.09 to 0.33 ± 0.10 and 0.17 ± 0.08 to 0.24 ± 0.08, respectively, whereas very low heritability estimates for both measures were detected at weaning; and (ii) alpha diversity at week 15 had strong and negative genetic correlations with BF, - 0.53 ± 0.23 to - 0.45 ± 0.25, as well as with ADG, - 0.53 ± 0.32 to - 0.53 ± 0.29. CONCLUSIONS: These results are important for efforts to genetically improve the domesticated pig because they suggest fecal microbiota diversity can be used as an indicator trait to improve traits that are expensive to measure.

Methods of tagSNP selection and other variables affecting imputation accuracy in swine.

BACKGROUND: Genotype imputation is a cost efficient alternative to use of high density genotypes for implementing genomic selection. The objective of this study was to investigate variables affecting imputation accuracy from low density tagSNP (average distance between tagSNP from 100kb to 1Mb) sets in swine, selected using LD information, physical location, or accuracy for genotype imputation. We compared results of imputation accuracy based on several sets of low density tagSNP of varying densities and selected using three different methods. In addition, we assessed the effect of varying size and composition of the reference panel of haplotypes used for imputation. RESULTS: TagSNP density of at least 1 tagSNP per 340kb (~7000 tagSNP) selected using pairwise LD information was necessary to achieve average imputation accuracy higher than 0.95. A commercial low density (9K) tagSNP set for swine was developed concurrent to this study and an average accuracy of imputation of 0.951 based on these tagSNP was estimated. Construction of a haplotype reference panel was most efficient when these haplotypes were obtained from randomly sampled individuals. Increasing the size of the original reference haplotype panel (128 haplotypes sampled from 32 sire/dam/offspring trios phased in a previous study) led to an overall increase in imputation accuracy (IA = 0.97 with 512 haplotypes), but was especially useful in increasing imputation accuracy of SNP with MAF below 0.1 and for SNP located in the chromosomal extremes (within 5% of chromosome end). CONCLUSION: The new commercially available 9K tagSNP set can be used to obtain imputed genotypes with high accuracy, even when imputation is based on a comparably small panel of reference haplotypes (128 haplotypes). Average imputation accuracy can be further increased by adding haplotypes to the reference panel. In addition, our results show that randomly sampling individuals to genotype for the construction of a reference haplotype panel is more cost efficient than specifically sampling older animals or trios with no observed loss in imputation accuracy. We expect that the use of imputed genotypes in swine breeding will yield highly accurate predictions of GEBV, based on the observed accuracy and reported results in dairy cattle, where genomic evaluation of some individuals is based on genotypes imputed with the same accuracy as our Yorkshire population.

Estimation of linkage disequilibrium in four US pig breeds.

BACKGROUND: The success of marker assisted selection depends on the amount of linkage disequilibrium (LD) across the genome. To implement marker assisted selection in the swine breeding industry, information about extent and degree of LD is essential. The objective of this study is to estimate LD in four US breeds of pigs (Duroc, Hampshire, Landrace, and Yorkshire) and subsequently calculate persistence of phase among them using a 60 k SNP panel. In addition, we report LD when using only a fraction of the available markers, to estimate persistence of LD over distance. RESULTS: Average r2 between adjacent SNP across all chromosomes was 0.36 for Landrace, 0.39 for Yorkshire, 0.44 for Hampshire and 0.46 for Duroc. For markers 1 Mb apart, r2 ranged from 0.15 for Landrace to 0.20 for Hampshire. Reducing the marker panel to 10% of its original density, average r2 ranged between 0.20 for Landrace to 0.25 for Duroc. We also estimated persistence of phase as a measure of prediction reliability of markers in one breed by those in another and found that markers less than 10 kb apart could be predicted with a maximal accuracy of 0.92 for Landrace with Yorkshire. CONCLUSIONS: Our estimates of LD, although in good agreement with previous reports, are more comprehensive and based on a larger panel of markers. Our estimates also confirmed earlier findings reporting higher LD in pigs than in American Holstein cattle, especially at increasing marker distances (> 1 Mb). High average LD (r2 > 0.4) between adjacent SNP found in this study is an important precursor for the implementation of marker assisted selection within a livestock species.Results of this study are relevant to the US purebred pig industry and critical for the design of programs of whole genome marker assisted evaluation and selection. In addition, results indicate that a more cost efficient implementation of marker assisted selection using low density panels with genotype imputation, would be feasible for these breeds.